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ABSTRACT 

A representative sample of unevolved early B-type stars in nearby OB associations and the field is analysed 
to unprecedented precision using NLTE techniques. The resulting chemical composition is found to be more 
metal-rich and much more homogeneous than indicated by previous work. A rms scatter of '--^10% in abun- 
dances is found for the six stars (and confirmed by six evolved stars), the same as reported for ISM gas-phase 
abundances. A cosmic abundance standard for the present-day solar neighbourhood is proposed, implying mass 
fractions for hydrogen, helium and metals of X = 0.715, F = 0.271 and Z = 0.014. Good agreement with solar 
photospheric abundances as reported from recent 3D radiative-hydrodynamical simulations of the solar atmo- 
sphere is obtained. As a first application we use the cosmic abundance standard as a proxy for the determination 
of the local ISM dust-phase composition, putting tight observational constraints on dust models. 
Subject headings: stars: abundances — stars: early-type — stars: fundamental parameters — ISM: abundances 
— dust, extinction — solar neighbourhood 
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1. EVITRODUCTION 

The Sun is unique among the stars because independent in- 
dicators allow its chemical composition to be constrained with 
a precision unmatched for any other star This can be done 
by spectroscopic analysis of its photosphere and by measure- 
ment of solar wind and solar energetic particles. Solar nebula 
abundances can be determined from CI chondrites, which are 
unaltered since the formation of the system. The wealth of 
information established the Sun as the principal standard for 
the chemical composition of cosmic matter (e.g. Grevesse 
& Sauval 1998, GS98; Holweger 2001; Asplund et al 2005, 
AGS05). However, is a 4.6 Gyr old star indeed representative 
of the cosmic matter in its neighbourhood^ at present? 

Ideal indicators for pristine abundances are unevolved early 
B-stars of spectral types B0-B2. Slowly rotating stars are pre- 
ferred as their photospheres should be essential ly unaffected 
by mixing of CN-processed material ( iMaeder & Mevne^ 
I2OOQ) . The atmospheres of early B-stars are also unaf- 
fected by atomic diffusion that gives rise to p ecuharities o f 
metal abundances in many later-type stars (e.g.'Smith"1996|). 
A major practical advantage is also their relatively simple 
photospheric physics, which is represented well by classical 
model atmospheres, unaffected by complications such as stel- 
lar winds or convection. 

As a consequence, early B-stars in the solar neighbour- 
hood were subje ct of several NLTE studi es in the past (e.g. 



yies & Lambert 1992; KiHan 1992, 1994; Cunha & Lambert 
1994; Daflon et aL .1999, ,2001a.b. 2003 ; Lyubimkov et al.. 
2004l2005h . Overall, they found a wide range of abundances, 
by about a factor ~10, and an average metallicity of only 
^2/3 solar (GS98). Hence, the impression arose that the solar 
neighbourhood is chemically highly heterogeneous, and the 
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Sun anomalously metal-rich compared to young stars. 

Both findings are problematic in terms of Galactic chem- 
ical evolution. Dispersal of stellar nucleosy nthesis products 
increases the metallicity overtime (e.g. lChiappini et alj2003h 
and hydrodynamic mixing t ends to homogen ize the interstel- 
lar medium (ISM) locally (lEdmundsl 1 1 975h . Characteristic 
timescales for homogenization are short, ranging from 10^- 
10** yrs on scales of 100-lOOOpc (Roy & Kunth 1995). 

In contrast to the young stars the interstellar gas shows a 
high degree of che mical homogeneity in the solar neighbour- 
hood (Sofi£j |2004l) . with the rms scatter of mean abundances 
being ^^10%. However, the ISM gas phase is not suitable 
as a tracer for cosmic abundances because of selective deple- 
tion of elements onto dust grains. Here we reinvestigate the 
conundrum of inhomogeneous stellar vs. homogeneous ISM 
gas-phase abundances in the solar neighbourhood, motivated 
by the finding of homogeneous B-star abundances for carbon 
(see detailed analysis by Nieva & Przybilla 2008, NP08). 

2. SAMPLE ANALYSIS 

Six bright and apparently slow-rotating early B stars in the 
solar neighbourhood - randomly distributed in OB associa- 
tions and in the field, and covering a wide range of stellar pa- 
rameters - were observed in early 2005 at ESO/La Silla, using 
FEROS on the 2.2 m telescope. Spectra with broad wave- 
length coverage and resolving power A/AAw48 000 were 
obtained, at very high-S/N (up to ^-^800 in the B-band). 

The quantitative analysis of the sample stars was car- 
ried out following the hy brid NLTE approach discussed by 
iNieva & Przvbiui (I20071 NP07) and NP08. In brief, line- 
blanket ed LTE mode l atmospheres were computed with AT- 
LAS9 (lKurucdll993h and NLTE line-formation calculations 
were performed using updated versions of our codes DETAIL 
and Surface. State-of-the-art model atoms were adopted 
(see Table [T]!, which allow atmospheric parameters and ele- 
mental abundances to be obtained with high accuracy. 

Multiple independent spectroscopic indicators were consid- 
ered simultaneously for the determination of the atmospheric 
parameters, effective temperature T^ff and surface gravity 
\ogg: all Stark-broadened Balmer lines and 4-6 ionization 
equilibria, of He I/II, C II/III/IV, O I/II, Ne I/II, Si II/III/IV and 
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TABLE 1 

Stellar parameters & elemental abundances 





HR6165 


HR 3055 


HR1861 


HR2928 


HR 3468 


HR5285 


Sp. Type 


B0.2 V 


BO III 


Bl IV 


Bl IV 


B 1.5 III 


B2 V 


Association 


Sco Cen 


Field 


Ul 1 Ud 1 D 


Field 


Field 




d(pc) 


152±20 


438it57 


450di59 


481 ±63 


319±41 


155di20 


Tiff (K) 


32000±300 


31200±300 


27000d=300 


26300±300 


22900±300 


20800di300 


logg(cgs) 


yi on 1 n AC 


3.y5±0.05 


y1 1 O 1 A AC 

4. 12ihO.U5 


y1 1 C 1 A AC 

4. 15±0.05 


1 ^;a_la AC 
3.60±0.05 


A TO 1 A AC 

4.22±U.05 


5 (kms"') 


5±1 


8±1 


3itl 


3±1 


5±1 


3±1 


V sin ( (km s ) 


4±4 


29±4 


12itl 


14±1 


11±2 


18±1 


C(kms"') 


4±4 


37±8 




20±2 


20±1 






10.99±0.05 (20) 


10.94±0.05 (16) 


10.99±0.05 (14) 


10.99±0.05 (14) 


10.99±0.05 (14) 


10.99±0.05 (13) 


e(C II)'' 


8.27±0.14 (13) 


8.35±0.08 (10) 


8.32±0.10 (19) 


8.28±0.08 (18) 


8.36±0.10 (17) 


8.32±0.08 (20) 


£(CIII)'' 


8.31±0.11 (17) 


8.30±0.05 (7) 


8.36±0.03 (11) 


8.27±0.02 (5) 


8.47±0.04 (2) 


8.42±0.06 (2) 


£(ClV)'' 


8.34 (2) 


8.45 (2) 










£(N 11)'' 


8.l0±U.li (/3) 


/. / /itO.Uo (23 ) 


1 . /D±u. uy (oi ) 


o.UUihU. iz (oi ) 


~} m_Ln in 

/.yz±u. iu (jo) 


"7 "T/S-Ln no ^ A'l\ 
1 . /OihU.Uo (4/) 


S(U 1) 






S 87-1-0 (X\ 


9. 9.'\A-C\ OS 
o.ojZtU.UJ \J) 


Q Q9_|_fi a:i 


s 70-1-0 ns (i\ 


£(011)"= 


8.77±0.08 (51) 


8.79±0.10(41) 


8.74±0.11 (52) 


8.74±0.09 (46) 


8.80±0.09 (40) 


8.71 ±0.05 (45) 


£(NeI)f 


8.12±0.05 (2) 




8.12±0.08 (9) 


8. 11 ±0.09 (9) 


8.05±0.09(10) 


8.07±0.07 (14) 


£(Ne II)f 


8.14±0.07 (16) 


8.07±0.07 (8) 


8.08±0.09 (14) 


8.03±0.12 (8) 


8.06±0.03 (2) 




£(Mg II)E 


7.62±0.03 (3) 


7.60±0.01 (2) 


7.58±0.10 (6) 


7.56±0.03 (3) 


7.51±0.10 (6) 


7.50±0.05 (4) 


£(Sill)'' 






7.47±0.17 (2) 


7.56±0.08 (2) 


7.51±0.10 (5) 


7.22±0.13 (6) 


£(Silll)*' 


7.50±0.08 (8) 


7.48±0.08 (6) 


7.46±0.I1 (9) 


7.52±0.11 (8) 


7.53±0.17 (7) 


7.29±0.05 (9) 


£(SilV)h 


7.50±0.04 (10) 


7.51±0.18 (5) 


7.50±0.08 (3) 


7.48±0.14 (2) 


7.50±0.04 (2) 




£(Fe II)' 










7.38 (1) 


7.38 (1) 


£(Fe my 


7.38±0.12 (17) 


7.49±0.12 (5) 


7.44±0.09 (33) 


7.48±0.10(30) 


7.42±0.12(36) 


7.40±0.09 (32) 



£(E1) = log(El/H) + 12, with rms uncertainties and number of analysed lines in parentheses. NLTE model atoms: H: IPrzvbT lla & Butler (200^; " IPrzvbilfl 
doQ^ : '' Nieva & Przybilla (2006, 2008); "^ Przybilla & Butler (2001); Przybilla et al. (2000); ■= Becker & ButieJ {T98f),updated; ' Morel & Butlei' HM^, gf- 
values of Froese Fischer & Tachiev 1 2004) forNel; ° Przybilla et al. (2001); '' Becker & Budet ( 1990), extended & updated; ' Becker ( 1998); i Morel et al. (200^^ 



Fe Will. Als o, the observed spectral e nergy distributions were 
reproduced dNieva & Przybillall2006h . The redundancy helps 
to avoid systematic errors. The microturbulent velocity was 
determined by demanding that abundances be independent 
of line equivalent widths. Elemental abundances e(El), rota- 
tional velocity v sin / and macroturbulence ( were determined 
from fits to individual line profiles. The results are summa- 
rized in Table[T] Stellar parameters and He and C abundances 
are identical with those derived by NP07/NP08, except for 
HR5285, where consideration of additional ionization equi- 
libria indicated small revisions, though agreement is obtained 
within the mutual uncertainties. Spectral types and spectro- 
scopic distances are also given in Table [T] agreeing well with 
HIPPARCOS parallaxes (HR6165, HR5285) or with the as- 
sociation distance (HR 1861). The positions of the stars in the 
Teff-log^-plane are indicated in Fig.[Tl where a comparison 
with evolution tracks is made. An overview of the location of 
the stars in the solar vicinity is also given there. 
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Fig . 1. — Comparison with stellar evolution tracks IMevnet & Maedeil 
1200 31). Observed and predicted N/C ratios (by mass) are indicated. The inset 
shows the location of the sample stars as projected on the Galactic plane, with 
the solar neighbourhood (as considered here) schematically outlined. Crosses 
denote positions of BA-SGs from the control sample (see text). 



The uncertainties in the atmospheric parameters were de- 
termined from the quality of the simultaneous fits to all di- 
agnostic indicators. Statistical uncertainties for abundances 
were obtained from the individual line data (rms values). Sys- 
tematic errors in the abundances due to uncertainties in atmo- 
spheric parameters, atomi c data and the qualit y of the spec- 
ti-um are ~0.1 dex (NP08. lPrzvbilla et alj[2006l) . i.e. about as 
large as the statistical errors. 

Our analysis of each individual star differs from standard 
studies in two main respects: l) practically all (unblended) 
lines of the ion spectra are analyzed instead of a few se- 
lected 'good' lines thus avoiding selection effects, and ll) all 
parameter indicators (in particular the ionization equilibria) 
are closely matched simultaneously, which has never been 
achieved before. As a result, practically the entire observed 
stellar spectrum is repr oduced closely by the spectral syn- 
thesis, see NP07, NP08. iPrzvbilla et alj (l2008h for examples. 
This is facilitated by the use of critically evaluated data in 
the model atom construction and a (time-consuming) iterative 
approach for a precise determination of the stellar parameters 
(NP08). Less accurate photometric Teff-estimates as adopted 
in most previous work are avoided, as Teff-uncertainties are 
often the most important sources of systematic error in the 
abundance derivation, next to ill-chosen atomic data (NP08) 
and logg-uncertainties (NP07). 

3. CHEMICAL HOMOGENEITY OF THE SOLAR VldNITY 

The status of previous NLTE abundance studies of early B- 
stars in the solar neighbourhood (covering clusters, associa- 
tions and field stars) is illustrated in Fig.|2] A wide range of 
abundance values is found for most elements, typically span- 
ning ~ 1 dex (for comparison, such a range is bridged by the 
cumulative effec t of ~13Gyrs of Galac tochemical evolution, 
see e.g. Fig. 2 of lChiappini et al.ll2003h . Moreover, the abun- 
dance distributions peak in most cases at sub-solar values, in 
particular when referring to the solar composition of GS98. 
Exceptions are He, where most previous studies find values 
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Fig. 2. — Comparison of chemical abundance studies ( NLTE) of B-stars in the solar neighbourhood. Red bars: present wor k; full line and thick full line (for the 
same stars as in our work): 'KiUarf * 1993, 1 19941); dotted: iCunha & Lamberll (1994l),ICunha et al.l 12009) for Ne; short-dashed: 'GieT&Lamber? ^19920, excluding 
bright giants and SGs; long-dashed: Daflon et al. 1 1999, 2001a b, 2003); dot-dashed: Morel & Buder (2008); triple-dot-dashed: Lvubimkov et al. (2004, 20050. 
Bin width is a/2 of the individual studies. Iron abundances in LTE in all previous work. HR5285 is Si-peculiar and therefore excluded from the silicon mean 
abundance. Solar abundances (O): GS98 & AGS05 (lower values). The panel for C is reproduced from NP08 for completeness. See the text for details. 



on average larger than solar, and Ne (about solar, GS98, from 
two very recent studies; see also Lanz et al. 2008 for the case 
of Ar). Several of thes e older B-star studies were combined 
bv lSnow & Wit3 (fT996h and Sofia & Meyer (2001, SMOl, see 
Table |2|l to derive a reference composition, inevitably result- 
ing in sub-solar average values and a large rms scatter. The 
former discrepancy has since been largely removed from a re- 
evaluation of solar abundances (AGS05). However, the sta- 
tus quo in terms of Galactochemical evolution can only be 
understood by invoking and fine-tuning extra processes such 
as infall/outflow of material and local retention of supernova 
products by large amounts. 

On the other hand, our sample of early B-stars implies a 
high degree of homogeneity for elemental abundances in the 
solar neighbourhood, with a scatter of ~10%, and absolute 
values of about solar (GS98 and/or AGS05, see Fig. |2] and 
Table |2]i. The only exception is N, which is most sensitive 
to mixin g of the atmospheric lay ers with CN-processed mate- 
rial (e.g. lMaeder & Meynelll2000l) . In this case the pristine N 
abundance may be indicated by the 3 objects with the lowest 
value, implying a pristine N/C ratio of 0.31 ±0.05 (by mass; 
error bar accounts for additional uncertainties). 

Although our sample is small, we regard it as representa- 
tive for the early B-star population in the solar neighbour- 
hood. The stars sample the relevant portion of the H-burning 
phase of the objects in the HRD in terms of T^s and \ogg 
(see Fig. [TJ. They also sample one hemisphere of the so- 
lar neighbourhood (inset of Fig. (TJ, half of them located in 
OB associations and the other half i n the field. All six stars 
were analyzed by lKiliarJ (ll992ilT994<) before, typically span- 
ning the entire abundance range in her sample of 21 stars (see 
Fig. |2]l. We therefore also find a chance selection of stars 
with similar chemical composition for our sample unlikely. 
The wide abundance ranges found in previous work reflect the 
lower accuracy of the analyses, while shifts of the abundance 
distributions relative to each other reflect systematics, with 



different temperature scales being the most important among 
these. Our results are supported further by a control sample 
of six BA-type supergiants (BA-SGs, Fig. [Hi, for which mean 
values of e(0) = 8.80±0.02 and e(Mg) = 7.55±0.07 were de- 
rived using the same analysis metho dology as applied here 
(iPrzybiUa et al.ll2006t lFirnsteinll2006l) . 

The finding of chemical homogeneity for our sample is 
in excellent accordance with results from the analysis of 
the ISM gas-phase in the solar neig hbourhood (Sofia 2004 
and references therein) and with theory reg arding the effi- 
ciency of hyd rodynamic mixing in the ISM ( Edmundslll97l 
iRoy & Kunthl[T9 95). Excellent agreement is also obtained 
with e lemental abundances in the Orion nebula (Esteban et al] 
|2004 E04, see Table |2]i, with the exception of C, which may 
be a consequence of the atomic data used in the Orion anal- 
ysis (see NP08 for the stellar case) plus overestimated dust 
corrections. 

In the following we briefly investigate the impact of this 
cosmic abundance standard on important topics of contem- 
porary astrophysics. 

4. THE COSMIC ABUNDANCE STANDARD: FIRST APPLICATIONS 

In general, excellent agreement of our B-star abundances 
with solar values from recent 3D radiative-hydrodynamical 
simulations of the solar atmosphere (AGS05) is obtained. The 
oxygen value falls between GS98 and AGS05 values and neon 
is compatible with GS98. This opens up new perspectives 
in the ongoing discussion on helioseismic constraints, chem- 
ical abundances and the solar interior model as reviewed by 
lBasu& Antial(l2008h . 

Our cosmic abundance standard also facilitates a precise 
determination of dust depletion in the local ISM for the pri- 
mary constituents. The amount of material incorporated into 
dust grains is determined by the difference between our B-star 
abundances and the ISM gas-phase abundances, see Table |2l 
Accordingly, a composition poor in carbon but rich in oxygen 



4 



Przybilla et al. 



TABLE 2 

Chemical composition of different object classes in the solar neighbourhood and of the Sun 





cosmic standard 


Orion 




Young 


ISM 


ISM 




Elem. 


B stars - tliis work 


gas+dust'' 


B stars'^ 


F&G stars'^ 


gas 


dusfi 


Sun=/f 


He 


10.98±0.02/ ■■• " 


10.988±0.003 










10.99±0.02 


C 


8.32±0.03/209±15 


8.52±0.02 


8.28±0.17 


8.55±0.10 


8.15±0.06!! 


68±26 


8.52±0.06/8.39±0.05 


N 


7.76±0.05/ 58± 7 


7.73±0.09 


7.81±0.21 




7.79±0.03'' 




7.92±0.06/7.78±0.06 


O 


8.76±0.03/575±41 


8.73±0.03 


8.54±0.16 


8.65±0.15 


8.59±0.01' 


186±42 


8.83±0.06/8.66±0.05 


Ne 


8.08±0.03/120± 9 


8.05±0.07 










8.08±0.06/7.84±0.06 


Mg 


7.56±0.05/ 36± 4 




7.36±0.13 


7.63±0.17 


6.17±0.02J 


34.8±4.4 


7.58±0.05/7.53±0.09 


Si 


7.50±0.02/ 32± 1 




7.27±0.20 


7.60±0.14 


6.35±0.05i 


29.6±2.2 


7.55±0.05/7.51±0.04 


Fe 


7.44±0.04/ 28± 3 




7.45±0.26 


7.45±0.12 


5.41±0.04i 


27.3±2.7 


7.50±0.05/7.45±0.05 



" in units of Iog(EI/H)+ 12 /atoms per lO' H nuclei - computed from average star abundances (mean values over all individual lines per element, equal weight 
per line); ''E04; '^SMOl; ''difference between the cosmic standard and ISM gas-phase abundances, in units of atoms per 10* H nuclei; '''^ GS98/AGS05, 
photospheric values; S lSofial 120041) ; '' IMever et al] (19971 corrected accordingly to Jensen et al. ( 2007); ' Cartledge et al. l 2004h; J lCarfledge et al.l (2006) 



and refractory elements is indicated. 

Such studies were undertaken previously, using e.g. abun- 
dances of the Sun, of B stars and of young F&G stars (e.g. 
ISnow& Wit"3ll996l SMOl, see Table IS as proxies for the 
determination of the dust-phase composition, however with 
mixed success. In particular, B stars were rejected as reli- 
able indicators as the derived abundances of material in dust 
at that time were too low to produce the observed interstellar 
extinction. Our study revives B stars as proxies of the ISM 
dust-phase composition, and even more so because of the ex- 
tremely low abundance scatter compared to all other standards 
considered so far, except for the Sun. 

The present results imply tight observational constraints on 
dust models in terms of carbon abundance. The observed 
properties of dust grains, as inferred from the interstellar ex- 
tinction law, have to be produced by a rather small amount 
of carbon, posing a challenge to most dust models (see e.g. 
ISnow & WittlfT995h . We can carry out a n imp ortant consis- 
tency check, following ICartledge et al.l (l2006h : the O pre- 
dicted to be incorporated in grains from the observed Mg, Si 
and Fe dust abundances and a rudimentary dust model agrees 
with the derived O dust abundance within the mutual (small) 
uncertainties. For the rudimentary dust model we assume sil- 
icates to be predominantly MgSiOs, with only a small frac- 
tion of Fe bound in silicates and only a small fraction being 
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of olivine-like composition. The remaining Mg and Fe frac- 
tion is considere d to be in oxi de form (MgO, FeO, Fe203, 
Fe304), see e.g. iDraind (l2003l) for a discussion of observa- 
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